Energy transfer in finite-size exciton-phonon systems : confinement-enhanced 

quantum decoherence 
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Based on the operatorial formulation of the perturbation theory, the exciton-phonon problem is 
revisited for investigating exciton-mediated energy flow in a finite-size lattice. Within this method, 
the exciton-phonon entanglement is taken into account through a dual dressing mechanism so that 
exciton and phonons are treated on an equal footing. In a marked contrast with what happens 
in an infinite lattice, it is shown that the dynamics of the exciton density is governed by several 
time scales. The density evolves coherently in the short-time limit whereas a relaxation mechanism 
occurs over intermediated time scales. Consequently, in the long-time limit, the density converges 
toward a nearly uniform distributed equilibrium distribution. Such a behavior results from quantum 
decoherence that originates in the fact that the phonons evolve differently depending on the path 
followed by the exciton to tunnel along the lattice. Although the relaxation rate increases with the 
temperature and with the coupling, it decreases with the lattice size, suggesting that the decoherence 
is inherent to the confinement. 

PACS numbers: 05.60.Gg,71.35.-y,71.38.-k,63.22.+m 



I. INTRODUCTION 

In molecular lattices, understanding how excitons 
carry energy from one region to another is a key step for 
explaining many phenomena Example among many 
are amide-I excitons in a-helices that favor the transduc- 
tion of the chemical energy into mechanical work 
vibrons in adsorbed nanostructures that enhance surface 
reactions or promote quantum information transfer [TlT - 
[l5j . and Frenkel excitons in light-harvesting complexes 
that convert solar energy into chemical energy |16l - [l9| . 
A molecular lattice exhibits regularly distributed atomic 
subunits along which the energy of an electronic transi- 
tion, or of a high-frequency vibrational mode, delocalizes 
owing to dipole-dipole interaction. This gives rise to a 
narrow-band exciton whose eigenstates are superimposi- 
tions of local states. When the dynamics is governed by 
the exciton Hamiltonian, such superimpositions are co- 
herent since a phase relation is kept between the local 
states. This favors a wavelike motion of the exciton that 
propagates coherently. 

Unfortunately, the exciton does not evolve freely but it 
interacts with the vibrations of the host medium usually 
responsible for dephasing. Therefore, the fundamental 
question arises whether the energy delocalizes coherently 
or incoherently (2fj| . For instance, recent experiments 
revealed that the coherent nature of the exciton may ex- 
plain the remarkable efficiency of light-harvesting com- 
plexes [2l|, H3| ■ Similarly, in a peptide helix, it has been 
shown that when a localized amide-I mode is selectively 
excited, the energy transport rate drastically increases 
[H[ suggesting that a coherent energy flow takes place 
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between amide-I modes |24j . 

From a phenomenological point of view, dephasing 
limited coherent motion was described within stochas- 
tic models [25W30I ] . The main idea is that the lattice vi- 
brations induce random fluctuations of the site energies 
that behave as independent stochastic variables. Con- 
sequently, the phase of each local state randomly fluc- 
tuates so that the coherent nature of the exciton grad- 
ually disappears. A wavelike motion takes place in the 
short-time limit whereas a diffusion like motion occurs 
in the long-time limit. However, it turns out that the 
previous scenario is not so simple when a microscopic 
description is used to mimic the influence of the fluctu- 
ating surrounding. To proceed, most theories are based 
on either the Frohlich model [31j or the Holstein model 
32] that provide a general description of an exciton cou- 
pled with acoustic o r op tic phonons, respectively (see for 
instance Refs. |33l442j |). Depending on the model pa- 
rameters, the exciton properties exhibit different facets 
ranging from quantum to classical, from weak coupling 
to strong coupling, from adiabatic to nonadiabatic and 
from large to small polarons [431 ] . 

In that context, a quite different result from that de- 
rived from stochastic approaches was obtained recently 
for a narrow-band exciton coupled with acoustic phonons 
in an infinite lattice. Indeed, within the nonadiabatic 
weak-coupling limit, that is, when the exciton moves 
slower than the phonons, it has been shown that the 
exciton propagates freely as if it was insensitive to the 
phonons [44]]. This feature was established by using a 
time-convolutionless generalized master equation (TCL- 
GME) for describing the evolution of the exciton reduced 
density matrix (RDM) (4f| [46J. In an infinite lattice, the 
phonons behave as a reservoir insensitive to the exciton 
so that the Born approximation applies 47]. The in- 
fluence of the phonons is encoded in a time-dependent 
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relaxation operator. Up to second order, it involves 
exciton-phonon coupling correlation functions whose dy- 
namics is governed by free phonons. Consequently, be- 
cause acoustic phonons exhibit spatial correlations over 
an infinite length scale, the relaxation operator describes 
a fast dephasing-rephasing mechanism which prevents 
the exciton diffusivity. Note that similar results have 
been obtained for a donor-acceptor pair coupled with a 
solvent that exhibits spatial correlations [48]. 

In the present paper, the exciton-phonon problem is 
revisited for describing energy redistribution in a finite- 
size lattice, a question that is rarely addressed. However, 
characterizing size effects is of fundamental importance 
because in biology and in nanoscience, relevant structures 
have reduced dimensionality. In a confined environment, 
we are faced with a major problem because the phonons 
no longer behave as a reservoir [49| . Indeed, as observed 
in the TCL-GME that governs the evolution of the exci- 
tonic coherences, that is, the off-diagonal RDM elements 
between the vacuum and the one-exciton states, the con- 
finement favors quantum recurrences and strong memory 
effects in the exciton-phonon coupling correlation func- 
tions 50]. Therefore, the relaxation operator is an al- 
most periodic function so that the TCL-GME reduces to 
a linear system of differential equations with almost pe- 
riodic coefficients. In that case, parametric resonances 
give rise to an exponential growth of the RDM indicat- 
ing that the TCL-GME method breaks down. The Born 
approximation fails to capture the dynamics because the 
correlations between the exciton and the phonons are no 
longer negligible during the propagation. 

To overcome these problems inherent in the non- 
Markovian dynamics, different strategies can be used, 
such as the correlated projector method [51], the 
time-dependent projection-operator approach [521 ] . the 
effective-mode representation [53j ] and the quantum 
jumps approach (54[. Here, we use the operatorial for- 
mulation of the perturbation theory (PT) [H[ recently 
introduced for studying the dynamics of the excitonic 
coherences [56l - [58| . In the nonadiabatic weak-coupling 
limit, this method is a powerful tool for describing the 
spectral properties of the exciton-phonon system over a 
broad energy scale. It is particularly suitable for char- 
acterizing the coherence dynamics and is more accurate 
than TCL-GME. Within PT, the system dynamics is gov- 
erned by an effective Hamiltonian that explicitly accounts 
for exciton-phonon entanglement. Exciton and phonons 
are thus treated on an equal footing to go beyond the 
Born approximation. 

The paper is organized as follows. In Sec. II, the 
exciton-phonon Hamiltonian is described. Then, the re- 
sults provided by PT are summarized to derive the ef- 
fective exciton-phonon Hamiltonian. Finally, the exciton 
RDM is defined as the key observable required for char- 
acterizing the energy transfer. Its expression is estab- 
lished within PT. In Sec. Ill, the simulation of the RDM 
dynamics is carried out numerically to extract the time 
evolution of the exciton density. The results are finally 



discussed in Sec. IV. 



II. THEORETICAL BACKGROUND 

A. Finite-size exciton-phonon system 

In a finite-size lattice containing N sites x — 1, N, 
the exciton- pho non Hamiltonian is defined as (see for in- 
stance Refs. (13,11113) 

H = H A + H B + AH. (1) 

In Eq.([T]), the first term is the exciton Hamiltonian Ha 
that acts in the one-exciton subspace £ a- It refers to 
the dynamics of N coupled two-level systems with Bohr 
frequency ujq. Note that a two- level system represents 
for instance the first two vibrational states of a localized 
high-frequency intramolecular vibration. The Hamilto- 
nian Ha is expressed in terms of the bare hopping con- 
stant $ between nearest neighbor sites as (note that the 
convention ft = 1 will be used throughout this paper) 

N N-l 

H A = J2 "o\x)(x\ + *(k + 1><*I + \x)(x + 1|), (2) 

x — l x—l 

where \x) is the first excited state of the xth. two- level 
system. Owing to the confinement, the exciton eigen- 
states are stationary waves with quantized wave vectors 
K k = kn/L, with k = 1, .., N and L = N + 1, as 

N [J 

\k) = J2\/j s HKkx)\x). (3) 

x — l 

The corresponding eigenenergies = ujo + 2$ cos(-Kfe) 
form a symmetric ladder of N discrete energy levels that 
belong to a band centered on luq and whose width is ap- 
proximately 4$. In the eigenbasis, the exciton Hamilto- 
nian is thus defined as : Ha — SitLi 

In Eq.([T]), the second term is the phonon Hamilto- 
nian Hb — Y2p=i ^p a p a p that acts in the Hilbert space 
8b- It describes the external motions of the lattice sites 
that behave as point masses M connected via force con- 
stants W. The phonons correspond to N stationary 
normal modes with quantized wave vectors q p = pir/L, 
with p = 1,..,N. The corresponding frequencies are 
Q p = fl c sin(q p /2) with fl c = y/AW/M. The dynam- 
ics is described using standard phonon operators a' p and 
a p and the eigenstates are the well-known number states 
denoted \n p ) = \n\, tin)- 

The last term in Eq.([l]) stands for the exciton-phonon 
interaction AH = J2 P =i M p (a p + a p ) that acts in Ea ® 
£b- It is defined in terms of the coupling matrix M p 
that measures the strength of the interaction between 
the exciton and the pth phonon mode. The interaction 
yields a stochastic modulation of each two-level system 
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Bohr frequency by the lattice vibrations. In the exciton 
eigenbasis, the matrix elements of M p are expressed as 



Mpkk' — Vp(3p,k-k' +8p,k'-k~Sp,k+k' —5p 



,2L-k-k' 



), (4) 



where rj p = [(Eb^ p / L)(l — (f^p/Hc) 2 )] 1 / 2 is expressed in 
terms of the small polaron binding energy Eb — X 2 /W 
that involves the coupling strength \. Equation Q shows 
that AH favors exciton scattering from state \k) to state 
\k') via the exchange of a phonon p. The allowed transi- 
tions are specified by the selection rules M p kk' that 
generalize the concept of momentum conservation in a 
finite-size lattice. 

The Hamiltonian H will be used for characterizing 
exciton-mediated energy redistribution. To proceed, we 
shall focus our attention to the nonadiabatic (4$ -C O c ) 
weak-coupling (Eb <C $) limit, a common situation 
for vibrational excitons in molecular lattices. In that 
case, AH being a small perturbation, the unperturbed 
states \k,n p ) refer to a free exciton accompanied by free 
phonons. Although AH favors transitions between un- 
perturbed states, the key point is that there is no res- 
onance between coupled unperturbed states so that PT 
can be used to solve H. 



B. Perturbation theory 

The operatorial formulation of PT is based on the in- 
troduction of a unitary transformation U that diagonal- 
izes the Hamiltonian H = U HU* in the unperturbed 
basis (5(1[57]]. It is written as U — exp(5), where S is an 
anti-Hermitian generator that is non-diagonal in the un- 
perturbed basis. Within PT, S is expanded as a Taylor 
series with respect to AH so that the diagonalization is 
achieved at a given order (see appendix A) . Up to second 
order, the transformed Hamiltonian is written as 



N 

# = E [(wfc + &Jfc)|*)(fc| 



fig ® 



k)(k\ , (5) 



k=l 



(k) 

where H B is the Hamiltonian that governs the phonon 
dynamics when the exciton lies in the state fc), as 

](Tl p + 5^l pk )ala p . (6) 



H 



(fe) 



N 

E( 



Equations (|5|) and ^ show that in a state |fc), the en- 
ergy of an exciton u> k = u>k + Sui k is renormalized due to 
its coupling with the phonons. Similarly, the frequency 
Cl p k — fl p + SQ p k of phonon p accompanied by an exciton 
in the state \k) is renormalized by an amount 5Q pk . The 
energy corrections are thus defined as 



5uj k 



N N 

EE 

p=i k'=i 



M 2 pkk , 



UJ k - UJ k , 



N 



$Qpk = 



- u k , ) 



k' = l 



(uj k - LU k >) 2 - £1 2 



(7) 



The operator H defines the effective exciton-phonon 
Hamiltonian. Being diagonal in the unperturbed basis, 
its eigenvalues e k , np = &k+J2p=i n p&pk stand for the sys- 
tem eigenenergies up to second order in AH. The corre- 
sponding eigenstates are defined as l^fc,,^} = U^\k,n p ). 
Consequently, U provides a new point of view in which 
H no longer describe independent excitations but refers 
to entangled exciton-phonon states. The entanglement 
results from a dual dressing effect because a state \k) 
defines an exciton dressed by a virtual phonon cloud 
whereas a number state \n p ) describes phonons clothed 
by virtual excitonic transitions. The behavior of the en- 
ergy corrections has been studied in great details (56l - [58| . 
Owing to the dressing, we have shown that the excitonic 
eigenenergies are redshifted and the energy correction of 
a stationary wave \k) scales as 



8uj k 



-E B (l-2/L) 



^7 C0S( T } - (8) 



By contrast, the phonon frequencies are either redshifted 
or blueshifted, depending on the state occupied by the 
exciton that accompanies the phonons. Such a mecha- 
nism is well-accounted by the following law 



SQ 



pk 



cos( — ) 



sm(— )cos ( 



-2L' 

Sp,k 



2L> 
Sp. L — k 



(9) 



To conclude, let us mention that for a given temper- 
ature and a fixed coupling strength, we have introduced 
a critical length L* w 0.m 2 /(E B k B T) (k B is the Boltz- 
mann constant) so that for L < L* , PT correctly de- 
scribes the exciton-phonon dynamics. By contrast, PT 
breaks down for L > L* owing to the occurrence of quasi- 
resonance between unperturbed states. 



C. Transport properties 

According to the standard theory of open quantum 
systems [45l |46| , a complete understanding of the exci- 
ton dynamics is obtained from the knowledge of the RDM 
<r(t) = TrB[exp(—iHt)p(0) exp(iHt)], where Tr B stands 
for a partial trace over the phonon degrees of freedom. 
The time evolution corresponds to an Heisenberg repre- 
sentation with respect to H and the initial conditions are 
specified by the exciton-phonon density matrix p(0). 

To define p(0), let us mention that without any pertur- 
bation the lattice is assumed to be in thermal equilibrium 
at temperature T. Assuming ujq ^> ksT, each two-level 
system is in its ground state. This is no longer the case 
for the phonons that form a thermal bath described by 
the Boltzmann distribution p B — exp(— /3H b )/Zb, where 
Zb is the phonon partition function (/3 = l/fceT). Con- 
sequently, to study exciton-mediated energy flow, one as- 
sumes that the lattice interacts with an external source 
which brings the system in a state out of equilibrium. 
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This source is resonantly coupled with the two-level sys- 
tems, only, and it does not affect the remaining degrees 
of freedom. Moreover, it is supposed to act during a very 
short time scale so that the exciton is prepared adiabat- 
ically without any change in the quantum state of the 
phonons. To simplify the discussion, we consider a situa- 
tion in which an exciton is created on the site xq at time 
t = so that the initial density matrix is p(0) — pa® Pb 
where pa = l^oX^ol- Note that for amide-I exciton in 
biopolymers, such an excitation may result from the en- 
ergy released by the hydrolysis of ATP [2] or from charge 
neutralization upon electron capture by a protonatcd he- 
lix 0. 

In that context, in the local basis, the exciton RDM is 
expressed as 

a(x,x',t) = Tr B [p B {x \e lHt \x'){x\e- lHt \x )}. (10) 

Its diagonal elements P(x,t) = a(x,x,t) yield the exci- 
ton density that measures the probability to observe the 
exciton on site x at time t given that it was created on 
site :ro a t time t = 0. The exciton density is the cen- 
tral object of the present study. It provides information 
about the way the excitonic energy flows along the lattice 
after its initial implementation. Moreover, its knowledge 
allows us to characterize the influence of the phonons on 
the process of energy redistribution. The off-diagonal el- 
ements of the RDM measure the ability of the exciton 
to develop or to maintain coherent superimpositions of 
local states. They provide additional information for de- 



termining whether the energy redistribution is coherent 
or incoherent. 

For describing the time evolution of the RDM, we do 
not derive a GME. Instead, we apply PT that directly 
provides an expression of the RDM. This procedure can 
be summarized as follows (see appendix B). We first in- 
troduce U and diagonalize H in Eq. (ITU1) . Second, we use 
the fact that H is the sum of independent contributions, 
each contribution describing the exciton-phonon system 
when the exciton occupies a specific state |fc) (Eq.([5])). 
Then, we define the following Heisenberg representation 



17* (t) 



Finally, we introduce the ef- 



fective phonon density matrix 



exp 



-PH B - it(H ( B k) - H B 



Z 



(kk') 



(*) 



where 



Z { f\t)=Tr B e-P HB - i « H z 



-H 



(<=')n 



(11) 



(12) 



Strictly speaking, p^ k ' (t) is not a density matrix be- 
cause it yields complex values for the phonon population. 
However, it is isomorphic to p B with the correspondence 
[3Q p — >• (3fl p + i(SQ p k — Sfl p k')t. It thus yields averages 
equivalent to thermal averages and facilitates the calcu- 
lations. In that context, after simple algebraic manipu- 
lations, the RDM is rewritten as 



J 



a(x, x , t) 



fe=l ft' 



z B 



e -i{u k -U) h r)tji r 



P^ fe ° (*) <^o [ (— *) ] fc'> <fc' | C/lor'Xa; | t/t | fc > <fc | C7 fe (— t) |^o)l ■ (13) 



r 



The final step, quite fastidious, consists in expanding U 
as a Taylor series with respect to AH. By bringing to- 
gether the various terms, one finally obtains the second 
order expression of a(x, x' , t) that will be used in the fol- 
lowing of the text. Note that we felt it was unnecessary 
to give the expression of the RDM that exhibits about 
thirty different contributions. 



called zero order approximation of the RDM written as 

N N r,{kk'). 



CO 



( x ,x',t) = 



'(*) 



-i(Cj k -uiy)t 



k=l k' = l 

U\IU\~J 



Z B 



x (x \k')(k'\x'){x\k){k\x ). 



(14) 



III. NUMERICAL RESULTS 



To conclude, let us mention that in the weak-coupling 
limit, the exciton-phonon coupling mainly induces a 
rcnormalization of the system energies without signif- 
icantly modifying the quantum states. The exciton- 
phonon eigenstates basically correspond to the unper- 
turbed states so that the transformation U behaves as 
the unit operator in Eq. (|13p . One thus obtains the so- 



In this section, Eq. (fT5|) is used for studying the ex- 
citonic energy flow in a finite-size lattice. Despite its 
general nature, the previous formalism will be applied 
for investigating energy redistribution in a lattice of H- 
bonded peptide units, a system for which the parameters 
are well-known [2l-[l0l| (Note that this formalism can be 
used for describing many situations in which a narrow- 
band exciton interacts with the phonons of the lattice 
such as Frenkel exciton in molecular aggregates [l] or vi- 



brational exciton in adsorbed nanostructures [llHlJ])- In 
such a lattice, peptide units linked by H bonds are reg- 
ularly distributed. Each site contains an amide-I mode 
(C=0 vibration) that gives rise to a vibrational exciton. 
This exciton interacts with the H bond vibrations that 
form a bath of acoustic phonons at temperature T. To 
study the amide-I dynamics, the following parameters are 
used: oj = 1660 cm" 1 , W = 15 Nm" 1 , M = 1.8 x 1CT 25 
kg, fl c = 96.86 cur 1 and $ = 7.8 cm" 1 . To avoid 
PT breakdown, the coupling strength will vary around 
X = 10 pN [lC§. Of course, the model displayed in Sec. 
II is too simple to accurately describe vibrational en- 
ergy flow in a real biopolymer whose dynamics exhibits 
a tremendous complexity due to the large number of de- 
grees of freedom. In particular, this model is unable to 
account for the finite amide-I lifetime because it con- 
serves the number of amide-I exciton. Nevertheless, it 
involves ingredients that play a key role in interpreting 
specific experiments such as pump-probe spectroscopy in 
a-helices [5| and Electron Capture Dissociation in finite- 
size polypeptides 8]. Consequently, we do not claim that 
the model is relevant to explain in details the vibrational 
dynamics in a protein. Nevertheless, its interest lies in 
the fact that it provides a simple approach to promote 
the idea that the confinement modifies the exciton energy 
redistribution. 



A. Time evolution of the exciton density 

As illustrated in Fig. 1 for N = 9, T = 300 K, x = 10 
pN and xq = 5, the evolution of the exciton density 
is governed by several time scales whose duration de- 
pends on the model parameters. In the short-time limit 
(t < 20 — 30 ps), Fig. la shows that P(x, t) flows almost 
coherently along the lattice. Note that a purely coherent 
motion is illustrated by the thick lines in Fig. la, as dis- 
cussed in the Sec. IV. The exciton behaves as a confined 
wave that experiences reflections on the lattice sides. In- 
deed, after the initial excitation, two wave packets propa- 
gate on each side of the central site xq- At t w 0.4 ps, the 
wave packets have left the excited region and more than 
50% of the initial energy occur on the neigboring sites 
xq±1. At t « 2.1 ps, the wave packets reach the lattice 
sides where the exciton density reaches 0.36. Then, the 
wave packets are reflected and they propagate back to 
the central site. Consequently, 71% of the initial energy 
recur on the central site at t ~ 7.1 ps. Such a mechanism 
takes place almost periodically so that quantum recur- 
rences occur. Therefore, P(xq, t) shows a series of peaks. 
The first peaks take place at t = 11.5 ps, t — 18.6 ps and 
t = 23.0 ps, the corresponding amplitudes being equal to 
0.81, 0.74 and 0.53. 

Over a longer time scale, the coherent nature of the 
exciton tends to disappear (Fig. lb). The recurrences in 
P(xo,t) become less and less pronounced. For instance, 
when the exciton moves coherently, P(xo-, t) exhibits two 
strong recurrences at t = 47.4 ps and t = 105.9 ps, the 
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FIG. 1: Time evolution of the exciton density P(x,t) for 
N = 9, T = 300 K, x = 10 pN and x = 5. (a) Short-time evo- 
lution and (b) time evolution over a longer time scale. Note 
that thick lines in Fig. la refer to a coherent energy trans- 
fer that arises when the relaxation induced by the phonons is 
neglected (see the Sec. IV). 



corresponding amplitudes being 0.98 and 0.97 (the sec- 
ond recurrence is not drawn). By contrast, when the 
exciton-phonon coupling is taken into account, P(xo,t) 
remains smaller than 0.46 in the neighborhood of the first 
recurrence and it does not exceed 0.27 in the neighbor- 
hood of the second recurrence. In fact, P(x,t) exhibits 
several spectral components. It thus develops a rather 
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FIG. 2: Time evolution of the exciton density P(xo,t) for 
N — 9 and xq = 5. (a) Influence of the exciton-phonon 
coupling strength for T = 300 K and (b) influence of the 
temperature for ^ = 10 pN. 



complex time evolution quite similar to that of an al- 
most periodic function, the almost period being the sig- 
nature of the recurrences. However, as time increases, 
specific spectral components disappear so that, after ap- 
proximately 100 ps, P(x,t) reduces to a damped sine 
function (Fig. lb). These features clearly characterize a 
relaxation mechanism that affects the energy redistribu- 
tion over an intermediate time scale. 

Finally, in the long-time limit, all the spectral compo- 
nents vanish. The density converges toward a station- 
ary distribution that specifies a new equilibrium density 
P eq (x). With the parameters used in Fig.l, P eq (x) is 
not uniformly distributed. At t — 500 ps, we obtain 
P(xo,t) ps 0.2 whereas P(x,t) ps 0.1, Vx ^ x n , suggest- 
ing that the exciton keeps the memory of its initial state 
despite its interaction with the phonons. Note that for 
X = 10 pN, we have verified that the zero order approx- 
imation of the RDM (Eq.(JT3J|) yields a quite good esti- 
mate of the time evolution of the density. 
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FIG. 3: Influence of the initial position of the exciton a;o on 
the equilibrium distribution P eq (x) for N — 9, \ — 10 pN and 
T = 300 K. 

The rate at which the exciton reaches the equilibrium 
depends on both the coupling strength and the tempera- 
ture. This feature is illustrated in Fig. 2 that displays the 
time evolution of the density on the excited site P(xo,t) 
for different parameter values. As shown in Fig. 2a, the 
coupling \ enhances the relaxation rate. The larger x is, 
the faster the P eq (xo) is reached. To estimate the relax- 
ation time tr, we assume that P(xq, t) scales as a rapidly 
varying signal dressed by a smooth envelope function that 
describes an exponential decay toward equilibrium. In 
doing so, one successively obtains tr ps 100 ps, tr ~ 40 
ps and tr ps 20 ps for x = 8 pN, x = Yl pN and x = 16 
pN, indicating that tr oc 1/x 2 - Similarly, as displayed 
in Fig. 2b, the temperature also enhances the relaxation 
rate. One thus obtains tr ps 190 ps, tr ps 90 ps and 
tr ps 60 ps for T = 100 pN, T = 200 pN and T = 300 
pN. Such a behavior suggests that the relaxation time 
is inversely proportional to T. Note that the zero order 
RDM (Eq.(|14jl) yields results almost identical to those 
displayed in Fig. 2. 

B. Equilibrium distribution 

The influence of the initial position of the exciton on 
the equilibrium distribution is illustrated in Fig. 3 for 
N = 9, x = 10 pN and T = 300 K. Although P eq {x) 
is invariant under mirror symmetry, that is, P eq (x) — 
P eq (L — x), two distinct situations occur. When xq — 
L/2, P eq {x) is partially localized on the excited site. It 
is almost uniformly distributed over the sites x ^ xq 
(P eq (x) ps 0.1 Vir ^ Xq) whereas it is approximately two 
times larger on the excited site (P eg (a;o) ps 0.17). In 
fact, when x nearly vanishes or at very low temperature, 
P eq {x) exhibits an universal behavior. It depends only on 
the lattice size and its scales as P eq (x) ps (1 + 8 XXo )/L. 
By contrast, when xq is not located at the center of the 
lattice, P eq (x) partially localizes over the two sites Xo 
and L — xq. It remains almost uniform over the sites 
x 7^ xq and x ^ L — xq (P eq (x) ps 0.1). However, it 
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FIG. 4: Parameter dependence of the equilibrium distribution 
for N = 9 and xo = 5. (a) Influence of the exciton-phonon 
coupling for T — 300 K and (b) influence of the temperature 
for x = 10 pN. 



becomes approximately 1.5 times larger on the sites Xq 
and L — xq. One thus obtains P eq {xo) = P eq (L — xq) — 
0.146, 0.142, 0.139 and 0.138 for x = 1, 2, 3 and 4, 
respectively. Note that when x nearly vanishes or at 
very low temperature, we have observed the following 
universal behavior: P eq {xo) — P eq {L — xq) rj 3/2L and 

Peq{x) Wl/i for X ^ Xq and X ^ L — Xq. 

The equilibrium distribution is also quite sensitive to 
the coupling and to the temperature, as shown in Fig. 
4 for N — 9 and xq — L/2. When \ increases (Fig. 
4a), P eq (a;o) decreases whereas P eq (x) increases Vx ^ xq. 
Whatever the x values, P eq (x) scales linearly with Eb- 
Therefore, when x ~ 16 pN {Eb ~ 0.86 cm -1 ), the equi- 
librium density becomes almost uniform in the neighbor- 
hood of the excited site. One thus obtains P eq {xo) ~ 
P eq {x ± 1) « 0.12. Then, when x ~ 18 pN {E B « 1.08 
cm -1 ) a hole occurs in the exciton density at the center 
of the lattice, that is, P eq {xo) < P eq {x) Vx ^ x . This 
hole continues to deepen when x increases again. Note 
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that, as discussed in Sec. II. B, PT is valid provided that 
X remains smaller than a critical value. With the param- 
eters used in Fig. 4a, this critical value is approximately 
X ~ 20 pN {Eb ~ 1-4 cm -1 ). Above this value, un- 
physical results were obtained, that is, P eq {xo) < when 
X > 26 pN. As shown in Fig. 4b, the influence of the 
temperature is quite similar to that of the coupling. In- 
deed, as T increases, P eq {xo) decreases whereas P eq {x) 



increases Vx ^ xq ■ Note that the density evolves linearly 
with the temperature. The key point is that over a broad 
temperature range, the temperature slightly affects the 
equilibrium distribution. For instance, P eq (xo) decreases 
from 0.2 to 0.17 when the temperature increases from 
to 300 K. Consequently, the exciton density does not 
become uniform in the neighborhood of the excited site 
and it does not exhibit a hole at the center of the lat- 
tice. Finally, a numerical fit of the data has revealed 
that P eg (a;o) ~ 2/L - 3r}{B)E B k B T/n 2 c , where 77(B) is 
a slowly varying function of the adiabaticity B = 2$/O c 
quite close to unity. 



C. Linear Entropy and relaxation rate 

Let us now focus our attention on the relaxation pro- 
cess. As shown in Figs. 1 and 2, it is not easy to define 
a relaxation time because the exciton density exhibits 
several spectral components that decay according to dif- 
ferent relaxation rates. To overcome this difficulty, a sin- 
gle relaxation time must be extract from a more " global 
observable" . To proceed, we consider the linear entropy 
defined as S{t) = 1 — TrA[cr{t) 2 }. It provides a mea- 
sure of the missing information about the state of the 
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exciton when compared with a situation in which the ex- 
citon evolves freely. It accounts for the arrow of time and 
it describes the relaxation mechanism as a consequence 
of an information transfer "from the exciton toward the 
phonon bath" in the form of cxciton-phonon correlations 

The behavior of the entropy is shown in Fig. 5 for 
N = 9, x = 5, T = 300 K and for two x values. Two 
kinds of calculations have been carried out by introduc- 
ing the zero order entropy So(t) (dashed lines), defined in 
terms of the RDM <ro(t) (Eq. lfTl)) ). and the second order 
entropy 52 (t) (solid lines), defined in terms of the RDM 
a(t) (Ea. (IT3l) '). The entropy So(t) is a monotonic function 
that increases with time. Initially equal to zero, it first 
scales as t 2 in the very short-time limit. Then, it increases 
as time increases by following an exponential function 
that rises to a maximum. In the long-time limit, it con- 
verges to a constant value Soc, quite close to unity, that 
depends only on L and xo- When the entangled nature 
of the exciton-phonon eigenstates is taken into account, 
the time evolution of the linear entropy slightly changes. 
As previously, £2 (t) increases from zero to reach a con- 
stant value. In the intermediate-time limit and in the 
long-time limit, it still behaves as an exponential function 
that rises to a maximum. However, two main differences 
occur. First, in the very short-time limit, S^i) rapidly 
increases over a few tenths of ps. Then, the entropy no 
longer defines a monotonic function. Instead, it behaves 
as an increasing function that supports high-frequency 
small-amplitude oscillations. 

Such a behavior allows us to introduce a phenomeno- 
logical measure of the relaxation rate. To proceed, we 
works with the zero order entropy whose evaluation is 
must easier. We thus assume that it scales as So(t) ss 
5oo(l — cxp(— 2Tnt)), the factor 2 in the exponential ac- 
counting for the fact that So(t) involves the square of 
the RDM. Of course, this expression works quite well in 
the intermediate-time limit and in the long-time limit, 
but it fails in reproducing the short-time behavior of the 
entropy So defined, the parameter Tr represents an ef- 
fective relaxation rate that provides information on the 
time tr = 1 /Tr needed to the exciton to reach the equi- 
librium. Note that, in the weak-coupling limit, it turns 
out that there is no much difference between the relax- 
ation rates provided by the two entropies. 

The behavior of the relaxation rate Tr is illustrated 
in Fig. 6 for xq = 1. As shown in Fig. 6a, Tr in- 
creases linearly with E B . For instance, at T — 300 K, Tr 
ranges from 0.068 cm -1 (tr — 77.85 ps) to 0.153 cm -1 
(tr = 34.59 ps) when \ varies from 10 pN (E B — 0.33 
cm" 1 ) to 15 pN (E B = 0.75 cm -1 ). Similarly, Fig. 6b 
reveals that Tr increases almost linearly with the tem- 
perature over a broad temperature range. For instance, 
for x = 12 pN (E B = 0.48 cm" 1 ), T R extends from 0.015 
cm" 1 (tr = 361.90 ps) to 0.098 cm" 1 (tr = 54.10 ps) 
when T varies from 50 K to 300 K. Note that we have 
also verified that Tr increases almost linearly with the 
adiabaticity B provided that the nonadiabatic limit is 





0.30 -1 




0.25 - 




0.20 - 


E 




0.15 - 


a 




u 


0.10 - 




0.05 - 




0.00 - 



T=100 K 
T=300 K 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

E B (cm -1 ) 





0.12 -, 




0.10 - 




0.08 - 


E 



0.06 - 






u 


0.04 - 




0.02 - 




0.00 - 




50 



1 1 1 1 r 

100 150 200 250 300 



T(K) 



E 

o 



« 0.10 




0.05 - 



FIG. 6: Behavior of the relaxation rate F_r. (a) Influence of 
the exciton-phonon coupling for N = 9, xq = 1 and for two 
values of the temperature, (b) Influence of the temperature 
for N = 9, xq = 1 and for two values of the coupling, (c) 
Influence of the lattice size for xo — 1, \ — 10 pN and T = 300 
K. 



9 



reached. Similarly, Tr slightly depends on io- A rather 
fast relaxation occurs for xq — L/2 whereas a quite slow 
relaxation takes place when xq = 2 or xq = L — 2. For 
instance, for N = 9, \ = 10 pN and T = 300 K, we 
obtained T R = 0.087 cm" 1 (t r = 60.50 ps) for x = L/2 
and T R = 0.065 cm" 1 (tr = 81.50 ps) for x = 2. Fi- 
nally, the most surprising effect arises from the depen- 
dence of the rate with respect to the lattice size L. In- 
deed, as illustrated in Fig. 6c, T R decreases with L. It 
approximately scales as oc 1/Lso that, for \ — 10 pN 
and T = 300 K, it decreases from 0.068 cm" 1 (tr = 77.85 
ps) to 0.039 cm -1 (tr = 135.97 ps) when L varies from 
10 to 20. This result suggests that the relaxation mecha- 
nism is an intrinsic property of the finite-size lattice. The 
relaxation is thus enhanced by the confinement so that 
the shorter the lattice is, the larger is the relaxation rate 
. To conclude, let us mention that the previous data were 
exploited to extract an analytical expression of the relax- 
ation rate. Therefore, in the nonadiabatic weak-coupling 
limit, it turns out that the rate scales as 



r, 



a(x ) 



BE B k B T 



(15) 



where a(xo) ss 6.5 ± 1 slightly depends on the initial 
position of the exciton. 



IV. DISCUSSION 

As shown in the previous section, our numerical results 
contrast sharply with what happens in an infinite lattice 
with translational invariance. Indeed, in this latter situa- 
tion, spatial correlations in the phonon bath prevent the 
exciton diffusivity. The exciton develops a wavelike mo- 
tion so that it propagates coherently along the lattice as if 
it was insensitive to the phonons. In a confined lattice, a 
fully different behavior takes place. Although the exciton 
delocalizes coherently in the short-time limit, its inter- 
action with the phonons favors a relaxation mechanism 
when time elapses. Consequently, the coherent nature 
of the exciton density gradually disappears. The density 
finally converges toward an almost uniform equilibrium 
distribution that slightly keeps the memory of the initial 
exciton position. Of course, in agreement with what one 
expects from a physical point of view, the relaxation rate 
is enhanced by the coupling strength and by the tem- 
perature. However, more surprisingly, we have observed 
that an increase of the lattice size softens the influence 
of the phonon bath resulting in a slowdown in the relax- 
ation process. In other words, the relaxation previously 
evidenced is an intrinsic property of a confined lattice. 



A. Approximate expression of the exciton density 

To understand the numerical results, let us discuss the 
physics which derives from the expression of the exciton 



density. According to Eq. (|T0|) and by expanding the par- 
tial trace in the phonon number state basis, the density 
is rewritten as 



P(x,t) = ^ (n p \p B \n p ) \(m p ,x\e~ 



lHt \x ,n p )\ 2 



(16) 



So defined, P(x,t) characterizes a generalized quantum 
probability. It measures the probability to observe the 
exciton-phonon system in the factorized state \x,m p ) at 
time t provided that it was in the state \xo,n p ) at time 
t = 0. Because we are concerned with the exciton dy- 
namics, a sum over all possible phonon states at time t is 
performed and an average over the initial phonon state 
is realized. 

To evaluate P(x,t), one can take advantage of the 
fact that in the nonadiabatic weak-coupling limit, the 
exciton-phonon coupling mainly induces a renormaliza- 
tion of the system energies without significantly modify- 
ing the quantum states. As a result, the Hamiltonian H 
in Eq.([TB|) can be replaced by the effective Hamiltonian 
H (Eq.©) so that the density becomes 



P(x,t) = 



N 

E< 

k=l 



" u " kt (x\k)(k\x a ) 



(17) 



where the symbol < ... > stands for an average over the 
initial phonon number state. In accordance with the laws 
of quantum mechanics, before the average, the exciton 
density reduces to the square modulus of the probability 
amplitude that the exciton tunnels from xq to x dur- 
ing time t, the phonons evolving freely from their initial 
state \n p ) . This amplitude is the sum over the elementary 
probability amplitudes associated to the different paths 
that the exciton can follow to tunnel. A given path de- 
fines a transition through a stationary wave \k) and it 
exhibits three contributions. First, it involves the weight 
of the localized states \x) and \x$) in the stationary wave 
\k). Then, it depends on a phase factor that accounts 
for the free evolution of the exciton. Finally, it involves 
a second phase factor that refers to the quantum evolu- 
tion of the phonons dressed by the exciton in the state 
\k). In that context, by performing the average over the 
phonon degrees of freedom, the exciton density is finally 
expressed as 



N N 



p{x,t) « Y^Ys F ( kk '^ e ~ i{c>h ~ Clh ' )t 



k=l k' = l 

U\liJ\i\ 



x {x \k'){k'\x'){x\k){k\x Q ), 



(18) 



where F(kk',t) stands for the so-called decoherence fac- 
tor defined as 



F(kk',t) = /e^IW^-^OA = 

\ / Zr 



Zb 



(19) 



At this step, it is straightforward to show that Eq.fTTS]) 
corresponds to the diagonal element of the zero order 
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RDM a (t) displayed in Eq.(Q3J}. Provided that k ^ k' , 
F(kk' , t) defines a decaying function that tends to zero in 
the long-time limit [13] • It represents the way the phonon 
bath destroys the coherence of an excitonic state defined 
as the superimposition between two stationary waves 
and Therefore F(kk',t) accounts for quantum deco- 
herence and it is expressed as 



the density is approximately expressed as 



A' 



F(kk',t) = U- 



1 



P =i 



-/3n p -i(5si pk -sn pk ,)t ■ 



(20) 



According to Eq. (|18p , the time evolution of the exciton 
density P(x, t) results from the quantum interferences be- 
tween the different paths that the exciton can follow to 
tunnel from xq to x. Owing to these interferences, P(x, t) 
is expressed as the sum of different spectral components, 
a specific component being associated to a phase factor 
involving an excitonic Bohr frequency. In addition, each 
component is modulated by a decoherence factor that 
results form the fact that the phonons evolve differently 
depending on the path followed by the exciton. More 
precisely, when the exciton follows two different paths, 
the phonons develop two different quantum evolutions so 
that an additional phase factor occurs. Strictly speak- 
ing, this phase factor is not responsible for decoherence. 
However, because the phonons are initially described by 
a statistical mixture of number states, an average is re- 
quired. This average yields a sum over phase factors that 
interfere with the other giving rise to quantum decoher- 
ence. 

As shown previously [13], the decoherence factor be- 
haves as a gaussian function whose time evolution can 
be extracted from the short-time limit of Eq.(|20p. One 
finally obtains 



F(kk',t) 



-iJ2 P n p (n pk -a pk ,)t -r 2 kk ,t 2 



(21) 



where "hp = [exp(/3f2 p ) — 1] _1 is the average number of 
phonons of the pth. mode. The so-called decoherence rate 
Tkk 1 j inversely proportional to the corresponding relax- 
ation time Tkk' = 1 /Tfcfc' , is defined as 



kk' 



\ 



1 N 



(22) 



p=l 



where An^ = n p (n p + 1) accounts for the thermal fluc- 
tuations of the phonon number. 



B. Coherent dynamics in the short-time limit 

By combining Eqs. (fT5|) and (|2"Tj) . we are able to for- 
mally explain the behavior of the exciton density as fol- 
lows. In the short-time limit, that is, over a time scale 
shorter than the smallest relaxation time, each decoher- 
ence factor reduces to an undamped phase factor so that 



P(x,t) 



N 

fc=i 



(x\k)(k\x ) 



(23) 



where Qk defines the effective energy of the exciton in 
state \k), as 



N 

LJ k — LO k + Suj k + ^2 i 
p=l 



pk • 



(24) 



As illustrated by the solid lines in Fig. la, Eq.(|23]) shows 
that the exciton moves freely as if it was insensitive to the 
phonon bath. A coherent energy transfer takes place, this 
transfer being characterized by quantum recurrences that 
result from the confinement of the lattice. Nevertheless, 
although they preserve the coherent nature of the exciton 
propagation, the phonons affect the dynamical parame- 
ters that govern this wavelike motion through the concept 
of effective energy. Indeed, as shown in Eq. (|24]l . it is as if 
the exciton was characterized by temperature dependent 
eigenvalues Uk that exhibit two corrections. The first cor- 
rection, encoded in the parameters Suik, accounts for the 
fact that during its propagation the exciton is dressed 
by a virtual cloud of phonons. The second correction, 
^ p fipSflpk, results form the phase of the decoherence 
factors. 

In that context, a moment's reflection will convince 
the reader that Qk reduces to the small polaron energy 
in a state obtained by using a finite temperature 
mean field theory (see for instance Ref. [36]). Within 
this theory, the exciton-phonon coupling is partially re- 
moved by performing the so-called Lang-Firsov transfor- 
mation. This transformation provides a new point of 
view in which the elementary excitation is no longer an 
exciton but a small polaron, that is, a composite particle 
that describes an exciton dressed by a lattice distortion. 
Then, a mean field approach is invoked to treat the re- 
maining polaron-phonon interaction, i.e. an average is 
realized over the phonon degrees of freedom according to 
the phonon density matrix ps- In doing so, it turns out 
that the polaron energy in a state |fe) is expressed as 



ipol) 



= Lj -e B + 2^e- s ^ cos(— ), 



(25) 



where es = Eb{1 — 2/L) generalizes the concept of small 
polaron binding energy in a confined lattice [58j and 
whereas s(T) is the well-known temperature dependent 
band-narrowing factor. The band-narrowing factor re- 
duces to s(0) = 8Eb/3ttQ c at zero temperature where it 
scales as s(T) = AEsksT /ft 2 , at high temperature. In 
that context, by inserting Eqs.flS} and © into Eq.flMJ), 
it is straightforward to show that the effective energy u>k 
reduces to the polaron energy w^. p °^ in the weak-coupling 
limit, that is, for s(T) <C 1. This result is quite interest- 
ing because it establishes the link between the present ap- 
proach, based on a perturbative treatment of the exciton- 
phonon entanglement, and the mean field procedure usu- 
ally introduced in the small polaron theory. Within PT, 



11 



the temperature dependence of the effective energy orig- 
inates in the modification of the quantum evolution of 
the phonons owing to exciton-phonon correlations. 



components of P{x, t) have disappeared owing to quan- 
tum decoherence. The exciton density converges toward 
a stationary distribution defined as 



C. Relaxation mechanism 

Over a time scale longer than the shortest relaxation 
time but shorter than the longest relaxation time, the de- 
caying behavior of the decoherence factors can no longer 
be disregarded. Therefore, a relaxation mechanism takes 
place resulting in the decay of the different spectral com- 
ponents of the exciton density. Each component dis- 
appears according to a specific decoherence rate. As 
shown in Eq. (|22p . the temperature dependence of T^k' 
is encoded in the fluctuations of the phonon numbers. 
Because An p approximately scales as fcsT/Jlp, Tkk' in- 
creases linearly with the temperature. Note that this 
dependence differs from the standard expression of the 
relaxation rate that characterizes the dephasing of an 
open system coupled with a reservoir of harmonic oscil- 
lators [l|. Indeed, in most situations the temperature 
dependence of the rate originates in its dependence with 
respect to the average phonon number n p . At high tem- 
perature, both approaches yield a similar temperature 
dependence since An p ss n p sa k B T '/f2 p . This is no longer 
the case at low temperature since An p /h p = exp(/?f2 p /2). 
Moreover, Tkk' is proportional to Er through its depen- 
dence with respect to the phonon energy corrections (see 
Eq.©). From these dependences, we thus recover that 
both the temperature and the coupling strength enhance 
the relaxation mechanism, as observed in Fig. 2. In that 
context, it is obvious that the duration of the relaxation 
mechanism will be approximately given by the longest 
relaxation time tr = Max[Tkk'], that is, the invert of 
the smallest decoherence rate Tr = Min[Tkk'}- From 
Eg. (|2"2"1) . it turns out that the smallest decoherence rate 
is ri 2 . Therefore, by inserting Eq.© into Eq. ([2"2"l) . one 
obtains a typical expression for the relaxation rate, as 



T R = 



4BE B k B T 



1 



27tt 4 
16X3' 



(26) 



Equation (1261) , quite similar to the expression of the rate 
extracted form numerical data (see Eq.JTSJ), clearly evi- 
dences the role of both the exciton-phonon coupling and 
the temperature. In addition, it shows that the relax- 
ation is an intrinsic property of a confined lattice, the 
rate Tr decreasing with the lattice size L. Consequently, 
the larger is the lattice size, the smaller is the relaxation 
rate. Such an effect originates in the size dependence of 
the energy correction of the phonons that scales as 1/L, 
as shown in Eq.©. 



D. Equilibrium distribution 

Finally, in the long-time limit, that is, over a time scale 
longer than the longest relaxation time, all the spectral 



P eq (x) 



N 

£ 



\(x\k)(k\x )\ 



(27) 



Equation (|2"7|) represents a time independent "classical" 
probability, the word classical meaning that there is no 
longer quantum interferences. It thus reduces to the sum 
of the probabilities associated to the realization of each 
path that the exciton can follow to tunnel from xq to x. 
Consequently, P eq (x) no longer depends on the temper- 
ature nor on the coupling strength. It involves only the 
lattice size and the initial position of the exciton. It is as 
if the exciton kept the memory of its initial location so 
that, for xq = L/2, the equilibrium distribution satisfies 



Pe q {x) 



2/L if X = Xq 
1/L if X ^ Xq. 



(28) 



Otherwise, for xq ^= L/2, the equilibrium distribution is 
defined as 



_ ■ 1.5/ L if x = xq and x = L — xq 
P eq {x) - \ x j L otherwise . 



(29) 



In a quite good agreement with the results displayed in 
Fig. 3, P eq {x) is partially localized on the sites xq and 
L — .To whereas it is almost uniformly distributed over the 
other sites. It is symmetric with respect to the center of 
the lattice. At equilibrium, the average position of the 
exciton is equal to x = L/2 and the standard deviation 
reduces to Ax 2 = (L 2 + 2)/12 + x (x Q /L - 1). 

In fact, such a behavior can be explained as follows. 
After straightforward calculations, it turns out the zero 
order RDM in the exciton eigenbasis is expressed as 

a(k,k',t) ^ F(kk' ,1)6-^"-^ (x \k'){k\x ). (30) 

Because the decoherence factor F(kk',t) decays with 
time provided that k ^ k' , the coherences of the RDM 
gradually disappear. By contrast, the populations are 
time independent and they remain unchanged. Conse- 
quently, the equilibrium corresponds to a statistical mix- 
ture of stationary waves \k). The corresponding proba- 
bilities, that is, \(xo\k)\ 2 , define the probabilities to ob- 
serve the exciton in states \k) given that it occupies the 
state |xo). From the expression of the stationary waves 
(Eq.©), one easily understands the characteristic of the 
equilibrium distribution. Note that within the previous 
scenario, we do not recover that P eq {x) depends on both 
the temperature and the coupling strength (see Fig. 4). 
In fact, the influence of these parameters results from 
the entangled nature of the exciton-phonon eigenstates, 
entanglement that is neglected within the zero order ap- 
proximation of the RDM. 
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E. Linear entropy 

To conclude, let us discuss the behavior of the linear 
entropy. Because the entropy measures the missing in- 
formation about the state of the exciton, its time evolu- 
tion can be interpreted as follows. At time t = 0, the 
exciton is in a pure state |xo) so that the entropy van- 
ishes. Then, when time elapses, quantum decoherence 
takes place. The coherences of the RDM gradually dis- 
appear so that the exciton is finally described by a statis- 
tical mixture of stationary waves. Consequently, this re- 
laxation mechanism is accompanied by an increase of the 
linear entropy (see Fig. 5). The zero order entropy con- 
verges toward a maximum value Soo = 1 — J2k K^ol^)! 4 
defined as 

Q _ / 1-1.5/L i£x ?L/2 . . 

0ao ~\l-2/L j£xo=L/2. 1 ' 

However, does the missing information correspond to 
an information transfer "from the exciton toward the 
phonon bath" in the form of exciton-phonon correlations? 
I don't think so. In fact, at zero order, the missing in- 
formation on the exciton state is due to our ignorance 
of the phonon states. Indeed, because the phonons are 
in thermal equilibrium, they are described by a statisti- 
cal mixture of number states. Therefore, the increase of 
the entropy originates in the quantum decoherence that 
results from the average over the phonon degrees of free- 
dom. The missing information on the initial phonon state 
is thus converted into missing information on the exci- 
ton state. Such a conversion is unidirectional so that 
the entropy behaves as a monotonic function. By con- 
trast, at second order, Fig. 5 has revealed that the linear 
entropy exhibits high-frequency small-amplitude oscilla- 
tions. These oscillations account for the fact that the 
process of information transfer between the exciton and 
the phonons is partially bidirectional. Basically, a certain 
amount of information is transferred from the exciton to 
the phonon bath. Then, owing to the entangled nature of 
the exciton-phonon eigenstates, a part of this information 
is restored to the exciton ... and so on. Nevertheless, ow- 
ing to quantum decoherence, this information exchange 
gradually disappears and the linear entropy overall in- 
creases as time increases. 



V. CONCLUSION 

In the present paper, the exciton-phonon problem was 
revisited for investigating exciton-mediated energy redis- 
tribution in a finite-size lattice. To proceed, the dynam- 
ics was described using the operatorial formulation of 
PT. Within this method, exciton-phonon entanglement 
is taken into account through a dual dressing effect, the 
exciton being dressed by a virtual phonon cloud whereas 
the phonons are clothed by virtual excitonic transitions. 
In the nonadiabatic weak-coupling limit, exciton and 
phonons were treated on an equal footing so that we were 



able to overcome the problem caused by the fact that the 
phonons no longer behave as a reservoir. 

In that context, special attention has been paid to de- 
scribing the exciton density that provides information 
about the way the excitonic energy flows along the lat- 
tice. In a marked contrast with what happens in an infi- 
nite lattice, we have shown that the dynamics of the den- 
sity is governed by several time scales. In the short-time 
limit, the density propagates coherently as if the exciton 
was insensitive to the phonons. Then, over an interme- 
diate time scale, the exciton-phonon interaction favors 
quantum decoherence so that the coherent nature of the 
density gradually disappears. Finally, in the long-time 
limit, the density converges toward an almost uniform 
equilibrium distribution. 

From a physical point of view, the exciton density mea- 
sures the probability for the exciton to tunnel between 
lattice sites, the phonons evolving freely. Therefore, its 
behavior results from the quantum interferences between 
the paths that the exciton can follow to tunnel. The key 
point is that when the exciton follows two different paths, 
the phonons develop two different quantum evolutions. 
The density thus depends on phase factors that reduce 
to time decaying functions after performing an average 
over the phonon degrees of freedom. Based on the char- 
acterization of the linear entropy, we have shown that 
the corresponding decay rate is enhanced by the exciton- 
phonon coupling strength and by the temperature. More 
surprisingly, our study revealed that an increase of the 
lattice size softens the influence of the phonons result- 
ing in a slowdown in the decoherence process. Quantum 
decoherence is thus an intrinsic property of a confined 
lattice. 

To conclude, we would like to draw the reader's at- 
tention to the fact that the present formalism must be 
restricted to the study of finite-size lattices. Indeed, we 
have shown that quantum decoherence tends to disap- 
pear as the lattice size increases. From these observa- 
tions, we could extrapolate to conclude that the exciton 
moves coherently in an infinite lattice. In doing so, we 
thus recover the results established in Ref.[44]. Unfortu- 
nately, this is a misleading appearance because the origin 
of the decoherence is twofold. First, it results from the 
modification of the phonon energies induced by the ex- 
citon. Consequently, the phonon quantum states evolve 
differently depending on the state occupied by the exci- 
ton. Dynamical phases arise in the exciton RDM so that 
pure dephasing occurs when the average over the phonon 
degrees of freedom is performed (see Eq.(|19p). Then, as 
illustrated by the presence of the unitary transformation 
U in Eq. (TT3l) . the decoherence also originates in the en- 
tangled nature of the exciton-phonon eigenstates. In a 
confined lattice, the entangled nature of the eigenstates is 
negligible so that the first origin of the decoherence dom- 
inates, as detailed in the present paper. By contrast, in 
an infinite lattice, the opposite situation occurs because 
the phonon energy corrections vanish (see Eq.@). For- 
tuitously, in the weak-coupling limit, it turns out that 
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the exciton-phonon entanglement does not favor an in- 
coherent (diffusion-like) motion of the exciton so that a 
coherent (wave-like) regime remains. But this is not al- 
ways the case. For instance, when one considers the dy- 
namics of the RDM elements that measure the coherence 
between the vacuum and the one-exciton states, we have 
shown that quantum decoherence takes place in finite lat- 
tices [571 ] . owing to exciton- induced phonon energy shift, 
and in infinite lattices [42| , owing to the entanglement of 
the exciton-phonon eigenstates. In that context, forth- 
coming works will be devoted to establish a quite general 
approach to treat the exciton-phonon system whatever 
the lattice size. In doing so, we aim at determining the 
nature of the decoherence process when the lattice size 
evolves continuously from a small value to infinity. 



Appendix B: General expression of the exciton 
RDM 



By inserting the unitary transformation U , the coher- 
ence Eq. flTUl) is rewritten as 



a(x,x',t) = (Bl) 



Since H is a sum of independent contributions (see 
Eq.©), the coherence is expressed as 



Appendix A: Generator of the unitary 
transformation 

As detailed previously [5?} , the generator of the unitary 
transformation up to second order in AH is given by the 
following equations 



N N 



[H A + H B ,Sx]=AH, 

{H A + H B ,S 2 } = ^[S u AH] nd , 



(Al) 



where nd means the non-diagonal part of an operator 
in the unperturbed basis. From the expression of the 
exciton-phonon Hamiltonian Eq.(fT]), one easily obtains 

51 = ^ Z p a} p - Z^a p , 

p 

5 2 = ^Eppta^al- E^apCLp, 



cr(x,x',i) = cxp [-i(uj k - uJk')t] 



k=l k'=l 



Trr 



i\U*\k) 



PB (x Q \U^\k')e m B \k'\U\x') 

(k) 



(k\U\x ) 



(B2) 



Because the operator exp(— iH B k ^t) defines a unitary evo- 
lution and because [H B ,H B k ^] = 0, the partial trace in 
Eq. (IB2[) is rewritten as 



Tr f 



PBe-^-^^ixole-^B^Uh+^^k') 



(k'lUlx^ixlU^^ikle-^B^Ue+^B^l^) 



+ D pp' a W - D lp' a l' a p +J2 c p- ( A2 ) 
pp' p 

In the exciton eigenbasis, the various operators that enter 
the definition of the generator are defined as 



where 



(k\Z p \k') 
(k\E P p,\k') 
(k\D pp ,\k') 

(k\C p \k') 



{k\M p \k') 

U) k - U) k , + O p 

(k\Bpp>\k') 

U! k — &k' + fl p + ilpi 

{k\B pp ,\k') 

Wfc — LUk' + O p — flpi 

(k\Af\k') 

UJk — UJk' 



B 



!'!•' c^ZpiM.pl\ 



A p = --(ZlMp + MpZp) 
B™6 P p> + Bp pl {\ - 5 P 



Bpp/ 



(A3) 



(A4) 



At this step, let first define U k (t) = e +iH B 't Ue -m B 't 
as the Hciscnberg representation of the unitary transfor- 
mation U with respect to the Hamiltonian H B k \ Then, 
since p B = exp(— j3 H B ) / 'Z B , one obtains 



with 



PB k '\i) = e -' 3 ^- 4 *(^ ) -^' , )/y lU ' 



^/Z {kk '\t) 



4 fcfc,) (t) = Tri 



-PH B -it{H ( g ) -H'g' ) ) 



(B3) 



Finally, inserting Eq. (|B3[) into the previous expression of 
the partial trace and combining the results with Eq. (|B2[) 
yield the exciton RDM Eq.^ 
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